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ABSTRACT 

Supernova remnants interacting with molecular clouds are potentially exciting systems in which 
to detect evidence of cosmic ray acceleration. Prominent 7-ray emission is produced via the decay 
of neutral pions when cosmic rays encounter the nearby dense clouds. In many of the supernova 
remnants coincident with 7-ray sources, the presence of OH(1720 MHz) masers is used to identify 
interaction with dense gas and to provide a kinematic distance to the system. In this paper we use 
statistical tests to demonstrate that there is a correlation between these masers and a class of GeV- to 
TeV-energy 7-ray sources coincident with interacting remnants. For pion decay the 7-ray luminosity 
provides a direct estimate of the local cosmic ray density. We find the cosmic ray density is enhanced 
by one to two orders of magnitude over the local solar value, comparable to X-ray-induced ionization 
in these remnants. The inferred ionization rates are sufficient to explain non-equilibrium chemistry in 
the post-shock gas, where high columns of hydroxyl are observed. 
Subject headings: cosmic rays — masers — supernova remnants — shock waves 



1. INTRODUCTION 

It has long been speculated that the vast energy of su- 
pernova remnants (SNRs) is tapped to produce the bulk 
of Galactic cosmic rays. Those SNRs interacting with 
molecular clouds are excellent targets to detect evidence 
of cosmic ray acceleration. The nearby dense gas acts as 
a target for hadronic cosmic rays, producing 7-rays via 
the decay of neutral pions (e.g., Drury et al. 1994). 

Interest in cosmic rays from interacting SNRs has 
been reinvigorated by recent detections of 7-ray emis- 
sion at GeV- to multi-TeV-energies from supernova rem- 
nants. The HESS galactic plane survey yielded 14 new 
TeV detections, 7 of w hich are coincident with SNRs 
(jAharonian et al.) 120061 ). Targeted observations are now 
capable of reaching sufficient sensitivities to resolve the 
morp hology of 7-ray em ission associated with the rem- 
nant (jAharonian et al.i r2008a,bl. In particular, a TeV- 
energy counterpart to IC 443 was observed t o be spatially 
offset from the GeV-ener gy 7-ray source (I Albert et al.l 
I2007t lAcciari et"alll2009f l. The spatial variations with 
energy may indicate the diffusion of cosmic rays accel- 
erated by the su pernova remnant o ut into the adjacent 
molecular cloud ([Torres et ani2008[) . 

Interacting SNRs represent a promising class of 7-ray 
sources which are likely to be uncovered in increasing 
numbers by the current generation of 7-ray observatories. 
However, to date the identification of 7-ray counterparts 
to SNRs has been inhibited by the poor spatial resolu- 
tion and sensitivity of 7 -ray telescopes (lEsposito et al.l 
Il996t [Torres et al.l [20031 ). Often, the association of 7- 
ray sources with interacting SNRs rests on supplemental 
evidence from common tracers of interaction with dense 
gas: broad molecular lines, bright at infrared to millime- 
ter wavelengths (e.g.. Reach et al. 2005), and the hy- 
droxyl(OH) maser at 1720 MHz (e.g.. Frail et al. 1994). 

OH masers, when detected only at 1720 MHz, are a 



particularly powerful diagnostic which is only associated 
with SNR shocks. While these "SNR masers" are only 
detected in ~10% of remnants, they unambiguously trace 
cool (25-150 K), dense (~10^ cm~'^) gas in the wake of 
non-dissociative shocks and per mit a direc t measure of 
magnetic field strength ( Lockett et al.|[l999l ) . Prominent 
interacting SNRs W28, W44 and IC 443 have been noted 
for both bright masers and coincident EGRET sources, 
leading to the suggest ion that masers may trace cosmic- 
ray acceleration sites (jClaussen et al.l[r997| ). 

In this letter we address the nature of interacting SNRs 
with 7-ray counterparts: First, we show that an excel- 
lent correlation exists between SNR masers and a sub- 
set of GeV- and TeV-energy sources toward SNRs. This 
strengthens the argument for a hadronic origin for the 
observed 7-rays. Second, if pion decay is indeed the 
dominant emission mechanism, the implied cosmic ray 
enhancement is sufficient to explain the increased ion- 
ization needed to produce h igh columns of OH in the 
post-shock gas ()Wardlelfl999il . This suggests there is a 
more intimate relationship between SNR masers and 7- 
ray sources than merely being complimentary tracers of 
dense gas interaction. 

2. 7-RAY EMISSION FROM INTERACTING SNRS 

First, we assess 7-ray sources which may be directly 
associated with supernova remnants (hereafter, "7-ray 
SNRs). To identify potential 7-ray counterparts to 
SNRs we have drawn from three catalogs: EGRET 
sourc es dTorres et a l. 2003), the Fermi Bright Source 
List ()Abdo et all 12009) and the TeVCat^. It is not 
trivial to determine the origins of each 7-ray source. 
Potential associations can be confused by the presence 
of pulsars or pulsar wind nebulae, which are also ca- 

^ TeVCat is an onUne catalog of very-high-energy (>50 GeV) 
7-ray sources. It can be accessed at http: / /tevcat. uchicago.edu . 
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TABLE 1 

Known SNR Masers and Coincident 7-hay Sources 



I 


b 


SNR 


Diameter 


Distance 


eV 


O-GeV 




"TeV 


Ref. 








(') 


(kpc) (ergs 


s-1) 




(ergs s~^) 






Group A 


6.4 


-0.1 


W28 


42 


2.0 


4.8e35 


2.1 


1.5e33 


2.7 


1,2 


34.7 


-0.4 


W44 


30 


2.5 


4.1e35 


1.9 






1 


49.2 


-0.7 


W51 C 


30 


6 


1.7e36 




1.5e33 




3,4 


189.1 


+3.0 


IC 443 


50 


1.5 


1.0e35 


2.0 


1.2e33 


3.1 


1,5,6 


Group B 


0.0 


+0.0 


SgrA East 


2.5 


8.5 


1.2e37 


1.7 


1.3e36 


2.2 




5.7 


-0.0 


9 


3.2 






0.8e33 


2.3 


2 


8.7 


-0.1 


W30 


45 


3.9 






1.7e34 


2.7 




337.8 


-0.1 


Kes 41 


5 


12.3 


4.2e36 


2.5 








348.5 


+0.1 


CTB 37A 


10 


11.3 


3.5e36 


2.3 


5.3e33 


2.3 


7 


359.1 


-0.5 




10 


5.0 


1.2e36 


2.2 


7.5e32 


2.7 


8 


Group G 


1.0 


-0.1 


Sgr D SNR 


8 


8.5 












1.4 


-0.1 




10 


8.5 












5.4 


-1.2 


Duck 


35 


5.2 












9.7 


-0.0 




11 


4.7 












16.7 


+0.1 




4 


2/14 












21.8 


-0.6 


Kes 69 


20 


5.2 












31.9 


-0.0 


3C 391 


8 


9 












32.8 


-0.1 


Kes 78 


20 


5.5/8.5 












337.0 


-0.1 


CTB 33 


3 


11 












346.6 


-0.2 




8 


11 












348.5 


-0.0 




10 


13.7 












349.7 


+0.2 




2 


>11 












357.7 


+0.3 


Square 


24 


6.4 












357.7 


-0.1 


Tornado 


5 


>6 












Note. 


— EGRET luminosity from 0.04-6 GeV and HESS or VERITAS luminosity from 0.3-10 TcV, 



wh ere a is the photon inde x. REFERENCES. - (1) [Esposito ct al. (1996), (2) Aharonian ct al. (20083), 
lAbdo et al.l (120091). (4) IFeinstein et al.l (120091). (5) lAcciari ct all I i2009l) . (6) lAlbcrt et al.l (i2007i) . (7) 



I Aharonian et a .ll^OOSaD . (8) lAharonian et a .1 (120061) 



pable of producing h igh-energy 7- r ay em ission. Some 
sources reported by iTorres et alj (|2003[ ) as potential 
EGRET counterparts to SNRs have since been associ- 
ated with other astrophysical sources. The source 3EG 
J1410-6147, coincident with SNR G312.4-0.4, has re- 
cently been associated with the young pulsar J1410- 
6132 with variability detected by AGILE (jO'Brien et all 
120081 ). The source 3EG J1824-151 4 has been con- 
firme d as the microquasar LS 5039 dAh aronian et al.l 
I2005D ruHng out a n association with SNR G16.8-1.1 
( Torres et al.l [20i51 . Sources 3EG J2016+3657 and 



3EG J2020+4017, coincident with SNRs G74.9+1.2 and 
G78. 2+ 2.1 respectively, ar e identified as extragalactic 
blazars (|Ivudin et al.ll2007f ). Finally, the EGRET coun- 
terparts for SN Rs G180.0-1.7, G 355.6+0.0, G39. 2-0.3 
and G74.8+1.2 (jTorres et al.' 20031 do not appear in the 
Fermi bright source catalog (Abdo ct al. 2009), and we 
therefore do not include them as detections in our anal- 
ysis. 

We include all remnants for which an association with 
the coincident 7-ray source has not been ruled out. There 
are 26 identified SNRs with 7-ray coincidences: 7 are 
young remnants (<1 kyr), 12 are SNRs with evidence 
of interaction with dense clouds, and 7 are unclassified 
remnants. Of the 12 identified SNRs which are interact- 
ing with dense gas, all but two (MSH 11-61A and W49B) 
have detected SNR masers. The presence of masers gives 
several advantages: (1) masers signpost interaction with 
dense (10^ cm~'^) clouds, which will enhance the pion 
decay signature, (2) the velocity of the maser gives a 
kinematic distance, and (3) an established velocity al- 
lows the adjacent cloud to be isolated from confusing 



Galactic emission along the line of sight. 

Table [1] lists all SNRs with masers and the properties 
of potentially associated 7-ray emission. In the first five 
columns the Galactic coordinates, name, diameter and 
kinematic distance arc listed for each remnant. The 7-ray 
luminosity for all derived from the reported ~100 MeV- 
10 GeV and >1 TeV fiuxes in columns 6 and 8. Spectral 
indices for are given in columns 7 and 9. References for 
detections are given in the last column. 

The ten SNRs with coincident 7-ray sources are di- 
vided based on the certainty of their association: Group 
A includes four SNRs for which 7-rays are established as 
related to the SNR; Group B includes six SNRs with co- 
incident 7-ray sources which have not been attributed to 
other astrophysical phenomenon (pulsars, blazars, etc.) 
but for which an association with the SNR is less than 
certain; Group C lists the sixteen SNRs with masers 
which do not yet have detected 7-ray sources. Given 
that both masers and 7-rays are detected for only 10% 
of SNRs, the large number of coincident detections makes 
an association between SNR masers and 7-ray emission 
as tracers of interaction quite plausible. 



2.1. Contingency Table Analysis 

To explore the correlation between SNR masers and 
7-ray sources, we use a contingency table analysis 
to test the null hypothesis that there is no associa- 
tion between the two groups. Here we include all 
remnants which have been searched for SNR masers 
(iFrail et all 119961: I Green et all 119971: iKoraleskv et al. 
119981: lYusef-Zadeh et all 119991: IHewitt k Yusef-Zadeh 
i2009i) and those 7-ray sources which have not been 
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TABLE 2 

Contingency Table: Presence of SNR Masers versus SNRs with 7-ray Sources 





7-ray SNRs 


No 7-ray Counterpart 








Expected 




Expected 




OH (1720 MHz) 




from 




from 




Masers 


Number 


Counts x'^ 


Number 


Counts x'^ 


Total 


Masers present 


10 


2.4 24.7 


14 


21.6 2.7 


24 


Masers absent 


14 


19.5 1.6 


185 


179.0 0.2 


199 


No observations 


2 


4.1 1.1 


41 


37.9 0.1 


43 


Total 


26 




240 




266 



Note. — We do not include those 43 SNRs which have not been observed for SNR iVIasers. 



clearly identified as having a non-SNR origin. 

Table [5] is the resulting contingency table. Three rows 
are given for SNRs with masers present, absent or not 
surveyed. Two columns divide SNRs with and without 
coincident 7-ray detections. Each of the two columns 
lists the "Number" of SNRs which meet both classifica- 
tions, the number "Expected from Counts" assuming the 
maser and 7-ray properties of remnants are completely 
independent, and "Contribution to x^" which gives the 
X^, (observed - expected) ^/(expected). For this test the 
total x^ is 31 with two degrees of freedom. This gives 
a probability of <0.0001% that SNR masers and 7-ray 
counterparts to SNRs are not correlated. 

We note that two of the expected cell frequencies are 
less than five. To be statistically rigorous we also apply 
the Fisher Exact Probability Test in addition to the x^ 
Test. The Fisher test considers all possible outcomes of 
the contingency table in order to determine the prob- 
ability of finding a correlation at least as strong as is 
observed. Here again we find a clear rejection of inde- 
pendence between 7-ray and maser detections in SNRs, 
with a probability of 0.002%. The two classes are clearly 
associated. 

A similar contingency table analysis was used to quan- 
tify an association between SNR masers and mixed- 
morphology remnants, a class of SNRs with shell-type 
radio morphologies with prominent thermal X-ra y emis- 
sion from their interiors (jYusef-Zadeh et ani2003[ ). Both 
classes of remnants are thought to result from interac- 
tion with dense gas and are strongly correlated. Of 
the identified 7-ray SNRs, 9 are classified as both 
mixed-morphology and maser-emitting, 2 are mixed- 
morphology without detected masers (SNRs W49B and 
MSH 11-61A), and only SNR G5. 7-0.0 has detected 
masers but no X-ray detection. Therefore, either clas- 
sification produces a nearly identically good correlation 
with 7-ray SNRs. Both X-rays and cosmic rays are ca- 
pable of providing the ionization needed to produce SNR 
masers. Detailed discussion is given in Section ??. 

An additional complication arises in that many SNRs 
have associated pulsars with their own wind nebulae that 
are also capable of accelerating particles to TeV energies. 
To account for this we separate the 7-ray SNRs into two 
categories based on whether a PWN is or is not also 
detected. This reduces the number of sources that fall 
under each classification, lowering the probability of in- 
dependence from the Fisher test to 1.2%, still significant 
enough to reject the null hypothesis. Furthermore, there 
is evidence that at least four SNRs (see Table [TJ " Group 
A") are associated with 7-ray emission from the SNR- 



cloud interaction, and not the PWN. Improved spatial 
resolution is needed to discriminate between 7-ray emis- 
sion from the PWN and SNR, but we note that there 
are no markedly different characteristics between the two 
groups of SNRs with and without PWN. 

2.2. Properties of ^ -ray SNRs with Masers 

For SNRs interacting with dense gas, 7-ray coun- 
terparts are likely to be dominated by emission from 
neutral pion decay originating from interactions be- 
tween hadronic cosmic rays and gas nuclei. In contrast 
to young SNRs, older interacting SNRs have little if 
any detected X-ray synchrotron emission. An inverse- 
Compto n scenario requires smal l magnetic field strengths 
<1 /iC ( Yamazaki et al.l l2006n whereas Zeeman split- 
ting of S NR masers meastires field strengths of order 
~1 mC (lYusef-Zadeh et"an [T996Hciaussen et al.llT997l: 
iBrogan et al.l I2000D . Electronic bremsstrahlung emis- 
sion has also been proposed, but the expectation would 
be f or 7-ray em ission to trace the radio shell morphol- 
ogy (iBvkov et a l. 2000). TeV emission from W28 and 
IC 443 shows no correlation with the radio shell, and 
an excellent correlation wi th dense gas (|Aharonian et al.l 
l2008bl : [Acciari et"aI1l2009f) . Forthcoming Fermi observa- 
tions will have sufficient resolution to resolve this ques- 
tion at GeV energies. 

The presence of a large reservoir of dense gas around 
the SNR, the expectation that cosmic rays are acceler- 
ated by SNRs, and the relatively older ages of this subset 
of interacting SNRs supports a pion decay origin for 7- 
ray counterparts. The morphology of 7-ray emission is 
seen to be well matched to that of the molecular cloud. 
In Figure [T] histograms of 7-ray luminosity show me- 
dians of 1.7x10^^ ergs s~^ and l.SxlO^'^ ergs s~^ for 
GeV and TeV detections, respectively. This is consistent 
with luminosity estimates if a significant fraction of the 
SN energy (10-20%) is diverted to cosmic-ra ys which are 
incid ent on the adjacent molecular cloud ( Drurv et al.l 
119941 ). The orders of magnitude differences in the lumi- 
nosity reflect the expected energy spectrum of acceler- 
ated particles. 

On average a hardening of the photon spectrum from 
GeV to TeV energies is also observed. The 7-ray spectral 
index, given in columns 7 and 9 of Table [U steepen at 
higher energies from a mean of aoeV ~ 2.1 to areV ^ 
2.6. It has been suggested for nearby SNR IC 443 that 
this spectral steepening results either from the pres- 
ence of bot h pion decay and a strong bremsstrahlung 
component (i Bvkov et al.ll200"ol ). or the diffusion of cos- 
mic rays at different energies through the dense cloud 
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Fig. 1. — Histograms of GeV (top) and TeV (bottom) luminosi- 
ties for SNRs with masers. Bin sizes of 1 and 0.5 are used. Char- 
acteristic luminosities of 1.7x10'^® ergs s~ 
are found, respectively. 



and 1.5x10'''' ergs s 



(|Torres et al.ll2008D . In addition to IC 443, the SNRs 
W28, Sgr A East, and G359.1-0.5 show this characteristic 
spectral hardening at higher energies, but CTB 37A does 
not. Future detailed spectral modeling of these sources 
can determine whether a spectral signature can be used 
to discriminate between different 7-ray emission mecha- 
nisms. 

3. ENHANCED IONIZATIONS VIA COSMIC RAYS 

The strong correlation observed between 7-ray-bright 
and maser SNRs may result from more than just 
both tracing interaction with high gas densities. The 
large post-shock OH abundances required for mas- 
ing to occur are not produced in chemical models of 
slow, dense shocks; instead all gas-phase oxygen is 
ra pidly converted to water b efore the gas cools to 50 
K (|Kaufman fc NeufeldlfTool . 

The large observed OH column densities are thought 
to be produ ced by dissoc iation of the abundant post- 
shock water (jWardldl 19991 ). In this model energetic elec- 
trons produced by ionizations in the molecular gas will 
excite H2, which will subsequently generate a weak flux 
of far-ultraviolet photon s which photo-dissociate water 
molecules producing OH (jPrasad fc Tarafdailll983[ ). The 
high flux of X-rays from the hot interior characteristic of 
mixed morphology SNRs has been shown to be a viable 
ionization mechanism (Yuscf-Zadch ct al. 2003). Cosmic 
rays and X-rays are notionally lumped together in this 
model. 

The 7-ray luminosity can be used to discriminate the 
relative importance of ionizations by cosmic rays. The 
7-ray spectrum resulting from pion decay of hadronic 
cosmic ray interactions in SNRs h as been described 
by many authors (Drux y ct al, 1994i ; Baring et al.lll99 9l: 
iBvkov et al.l[2000D . Here we follow 



AharonianI (| 199 If ). 



formulating the 7-ray flux as a function of the cosmic 
ray enhancement in the cloud adjacent to SNR, neglect- 
ing any possible gradients through the cloud as a simpli- 
fication. The 7-ray flux above energy E.y is given by the 



equation: 

F(> E,) ^ 3 X 10-" (^1^ 




ph cm s 



E- 



(1) 
is the 



where dkpc is the distance to the SNR in kpc, 
threshold above which the 7-ray flux is totaled, ujcb. is 
the ratio of cosmic ray density at the SNR to the local 
solar value assuming the same 7-ray emissivity, and M5 is 
the total mass of the SNR shell and adjace nt molecular 
cloud in units of 10^ Mq. We note that iDrurv et al.l 
()1994f ) established that spectral variations between 2.0 
and 2.7 changes this estimate by less than 20% . 

Assuming the observed 7-ray flux from the SNR is due 
to hadronic interactions with the surrounding dense gas, 
it is straightforward to roughly estimate the local energy 
density of cosmic rays, provided the distance and molec- 
ular environment are known. For EGRET and Fermi de- 
tections, E.y > 100 MeV, equation ([T]) can be rewritten 
to express the local cosmic ray enhancement ivcr- 



F^{> 100 MeV) 
7 X 10^'' ph cm^2 g- 



A/5 



(2) 



The cosmic ray ionization scales as the cosmic ray den- 
sity, so we obtain the ionization rate due to cosmic rays 
Ccfl, — ^CR C01 where ^0 i s the local cosm ic ray ioniza- 
tion rate of -4x10-1"^ s'^ (|WebbeHll998f ). 

In Table [3] we list the observed properties of SNRs 
from which the local enhancement of cosmic rays is 
estimated. We find cosmic ray densities are typically 
increased by 30-150 times that observed for quiescent 
molecular clouds. This provides an ionization rate suf- 
ficient to produce OH abundances of 10-^-10^^ in 
the post-shock gas. For comparison we also list the 
ionization rate from inte rior X-ray emission derived by 
lYusef-Zadeh et al.l (|200l . The two sources of ionization 
are generally found to be comparable. Either could be 
the dominant ionization mechanism depending on the lo- 
cation of the gas with respect to the interior X-ray emit- 
ting plasma and cosmic ray acceleration sites, and may 
vary on a source- by-source basis. 

We caution that our estimates of cosmic ray ionization 
rate are somewhat crude. Even if the dominant emission 
mechanism is hadronic, there may still be a significant 
bremsstrahlung component. Detailed spectral modeling 
with data from the Fermi 7-ray Space Telescope will 
permit a much better estimate of the local cosmic ray 
density. Observations of direct chemical trace rs of cosmic 
ray ionization r ates, such as Hj" an d He"*" (jOkal 120061 : 
iDalgarnol [200l llndriolo et ani2009f) . could be used to 
confirm the cosmic ray enhancements in the immediate 
environment of these interacting SNRs. 

4. SUMMARY & CONCLUSIONS 

We have explored the class of 7-ray sources which can 
be explained by enhanced cosmic ray densities at the 
interaction sites between SNRs and molecular clouds. 
By correlating SNR masers with known GeV and TeV 
sources, we have identified an emerging class of 7-ray- 
bright interacting SNRs. Of the 24 known Maser SNRs 
currently identified there are ten with GeV to TeV- 
energy 7-ray associations, and six with both. If 7-ray 
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TABLE 3 

Comparison of Cosmic Ray and X-ray Ionizations 



SNR Distance M^^ud F^(>100 MeV) Ccfl Cxray Rcf. 

(kpc) (105 Mq) (10-** cm-2 s-i) (10^16 s^i) (W^^ s-i) 

W28 2.0 0.5 74.2 3.4 2.3 1 

W44 2.5 3.0 88.9 1.1 4.5 2 

W51 C 6.0 1.9 40.9 4.4 8.8 3 

IC 443 1.5 0.1 51.4 5.5 3.6 4 



Not e. — References: (1) 
I I2009I 1. (4) IDickman et al. 



l lAbdo et al.ll200g|) 



Aharonian et al.l ll2008bf) . ('2') ISeta et al.l lll99gD . OI IFeinstein et al.l 
I I1992I) . Fluxes are taken from the Fermi Briglit Source List 



emission from these sources is largely due to hadronic 
cosmic rays, the enhanced local cosmic ray ionization 
rates in these clouds can explain the productio n of OH 
molec ules behind a C-type shock, suggested by IWardld 
(|1999f ). Furthermore, cosmic ray ionization is typically 
comparable to or dominant over ionizations from inte- 



rior thermal X-rays, though without detailed knowledge 
of the cosmic ray spectrum these results have large uncer- 
tainties. Interacting SNRs represent a promising class of 
7-ray sources which are likely to be uncovered by future 
7-ray observatories. 
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